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Abstract. We review how RHIC is expected to deepen our understanding of the spin structure of 
longitudinally polarized nucleons. After briefly outlining the current status of spin-dependent par- 
ton densities and pointing out open questions, we focus on theoretical calculations and predictions 
relevant for the RHIC spin program. Estimates of the expected statistical accuracy for such mea¬ 
surements are presented, taking into account the acceptance of the RHIC detectors. 


LESSONS FROM POLARIZED DEEP INELASTIC SCATTERING 


Before reviewing the prospeets for spin physios at the BNL-RHIC we briefly turn to 


longitudinally polarized deep-inelastio soattering (DIS) and what we have learned from 
more than twenty years of beautiful data [jT]]. To next-to-leading order (NLO) in the 
strong eoupling ct^, the DIS strueture funetion gj, whieh parametrizes our ignoranee 
about the nuoleon spin strueture, oan be expressed as 
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The symbol 0 denotes the usual eonvolution, and ACq,g are the perturbatively ealeulable 
eoeffioient funotions, whieh are known even up to next-to-next-to-leading order [@]. The 
A/, / = (q, q, g), are the spin-dependent parton distributions, defined as 

Af{x,ld)=f"^{x,il)-f~{x,n) , (2) 

where /+ (/^) is the number density of a parton type / with helioity “+” in a 
proton with positive helioity, oarrying a fraetion x of the proton’s momentum. Onoe they 
are known at some initial soale /Iq, their soale /i dependenoe is governed by a set of 
evolution equations 
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So far, the spin-dependent f i splitting funetions entering these evolution equations 
have been ealeulated up to NLO aoouraey [0 . 

Figure [H eompares the available information on gj (x, Q^) for DIS off a proton target to 
results of a typieal NLO QCD fit. From sueh types of analyses a pretty good knowledge 





FIGURE 1. Available information on gy{x,Q^) as collected by fixed-target experiments [|^ compared 
to results of a typical NLO QCD fit (solid lines). The indicated rectangular and triangular regions contain 
data which would not pass kinematical cuts of > 4GeV^ and W'^ > lOGeV^, respectively, usually 
imposed in all fits to unpolarized DIS data. 


of certain combinations of different quark flavors has emerged, and it became clear 
that quarks contribute only a small fraction to the proton’s spin. However, there is 
still considerable lack of knowledge regarding the polarized gluon density Ag, which is 
basically unconstrained by present data, the separation of quark and antiquark densities 
and of different flavors, and the orbital angular momentum of quarks and gluons inside a 
nucleon. With the exception of orbital angular momentum RHIC can address all of these 
questions as will be demonstrated in the following [Q]. 

There is also an important difficulty when analyzing polarized DIS data in terms 
of spin-dependent parton densities: compared to the unpolarized case the presently 
available kinematical coverage in x and ^ and the statistical precision of polarized DIS 
data are much more limited [[p. As a consequence, one is forced to include data into the 
fits from (x, 2^)-regions where corresponding fits of unpolarized leading-twist parton 
densities start to break down, see Fig. |T]. Data from RHIC, taken at “resolution” scales 
jl where perturbative QCD and the leading-twist approximation are supposed to work, 
can shed light on the possible size of unwanted higher-twist contributions in presently 
available sets of polarized parton distributions. 


















SPIN PHYSICS AT RHIC 


Prerequisites 


The QCD-improved parton model has been successfully applied to many high en¬ 
ergy scattering processes. The predictive power of perturbative QCD follows from the 
universality of the parton distribution and fragmentation functions which is based on the 
factorization theorem. To be specific, let us consider the inclusive production of a hadron 
H, e.g., a pion, in longitudinally polarized pp collisions at a c.m.s. energy \/S. The cross 
section can be written in a factorized form as a convolution of perturbatively calculable 
partonic cross sections describing the hard scattering ab cX and appropriate 

combinations of parton densities A/^ ^ and fragmentation functions embodying the 
non-perturbative physics: 

abc'' 


Here, F stands for any appropriate set of kinematical variables like the transverse mo¬ 
mentum pj^ and/or rapidity y of the observed hadron. The are the parton-to-hadron 
fragmentation functions. Their scale /i-dependence is governed by a set of equations 
very similar to (|]). The factorization scale jl, introduced on the r.h.s. of (^), separates 
long- and short-distance phenomena, ji is completely arbitrary but usually chosen to be 
of the order of the scale characterizing the hard interaction, for instance pj-. Since the 
l.h.s. of (Q) has to be independent of p. (and other theoretical conventions), any resid¬ 
ual dependence of the r.h.s. on the actual choice of p gives an indication of how well 
the theoretical calculation is under control and can be trusted. In particular, leading or¬ 
der (LO) estimates suffer from a strong, uncontrollable scale dependence and hence 
are not sufficient for comparing theory with data. Figure || shows a typical factoriza¬ 
tion scale dependence for various “high-py,” processes and experiments as a function 
of pj.. Clearly, the situation is only acceptable at collider experiments where one can 
easily reach py values in excess of 5 GeV. p^ values of the order of 1-2 GeV, accessi¬ 
ble at fixed-target experiments, are not sufficient to provide a large enough scale p for 
performing perturbative calculations reliably. For simplicity we have not distinguished 
between renormalization and initial/final-state factorization scales in which can be 
chosen differently. 

In practice, spin experiments do not measure the polarized cross section dAo/dT 
itself, but the longitudinal spin asymmetry, which is given by the ratio of the polarized 
and unpolarized cross sections, i.e., for our example above, Eq. (g), it reads 


yy _ dAo^ /dr 
do^/dT 


( 5 ) 


The unpolarized cross section do^/dT is given by Eq. (Q) with all polarized quantities 
replaced by their unpolarized counterparts. At RHIC one can also study doubly trans¬ 
verse spin asymmetries 0 but here we will focus on longitudinal polarization only. 
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FIGURE 2. Typical factorization scale dependence for various “high-pj.” processes and experiments as 
a function of pj. Shown is the cross section ratio for two choices of scale, pj and Pj/2. 


Accessing 

The main thrust of the RHIC spin program [Q] is to pin down the so far elusive gluon 
helicity distributions Ag{x, n). The strength of RHIC is the possibility to probe Ag(x, ^u) 
in a variety of hard processes 0], in each case at sufficiently large pj where perturbative 
QCD is expected to work. This not only allows to determine the x-shape of Ag(x, /i) for 
X >0.01 but also verifies the universality property of polarized parton densities for the 
first time. In the following we review the status of theoretical calculations for processes 
sensitive to Ag, experimental aspects can be found, e.g., in [^. 

The “classical” tool for determining the gluon density is high-p^ prompt photon pro¬ 
duction due to the dominance of the LO Compton process, qg ^ Jq- Exploiting this 
feature, both RHIC experiments, PHENIX and STAR, intend to use this process for 
a measurement of Ag. Apart from “direct” mechanisms like qg yq, the photon can 
also be produced by a parton, scattered or created in a hard QCD reaction, which frag¬ 
ments into the photon. Such a contribution naturally arises in a QCD calculation from 
the necessity of factorizing final-state collinear singularities into a photon fragmenta¬ 
tion function Dj. However, since photons produced through fragmentation are always 
accompanied by hadronic debris, an “isolation cut” imposed on the photon signal in 
experiment, e.g., a “cone”, strongly reduces such contributions to the cross section. 

The NEO QCD corrections to the direct (non-fragmentation) processes have been 
calculated in [|1] and lead to a much reduced factorization scale dependence as compared 
to EO estimates. In addition, Monte Carlo codes have been developed |^], which 
allow to include various isolation criteria and to study also photon-plus-jet observables. 
The latter are relevant for Ag measurements planned at STAR Bi- Since present 
comparisons between experiment and theory are not fully satisfactory in the unpolarized 








FIGURE 3. for prompt photon production in NLO QCD as a function of pj for different sets of 
parton densities. The “error bars” indicate the expected statistical accuracy 5A^^, Eq. (^, for the PHENIX 
experiment. Figure taken from [|^. 


case, in particular in the fixed-target regime, considerable efforts have been made to push 
calculations beyond the NLO of QCD by including resummations of large logarithms 
[^. It is hence not unlikely that a better understanding of prompt photon production 
can be achieved soon. Figure ^ shows as predicted by a NLO QCD calculation 
as a function of the photon’s transverse momentum py. The applied rapidity cut 
|7] I < 0.35 matches the acceptance of the PHENIX detector. The important result is that 
the expected statistical errors are considerably smaller than the changes in 
due to different spin-dependent gluon densities over a wide range of py. Hence RHIC 
should be able to probe Ag in prompt photon production. 5App may be estimated by the 
formula 


^4 = ^ 


( 6 ) 


where P is the polarization of one beam, J/f the integrated luminosity of the pp colli¬ 
sions, and the unpolarized cross section integrated over the p^-bin for which the 
error is to be determined. Unless stated otherwise, P = 0.7 and = 320 (800) pb^^ is 
used in Eq. @ for pp collisions at ^/S = 200 (500) GeV [|p. 

Jets are another key-process to pin down Ag at RHIC: they are copiously produced 
at Vs = 500 GeV, even at high pj, 15 < p^ < 50GeV, and gluon-induced gg and qg 
processes are expected to dominate in accessible kinematical regimes. Due to limitations 
in the angular coverage, jet studies will be performed by STAR only. As jet surrogates, 
PHENIX can look for high-py leading hadrons, such as pions, whose production pro¬ 
ceeds through the same partonic subprocesses as jet production. Hadrons have the ad¬ 
vantage that they can be studied also at \fS = 200 GeV and down to lower values in pj 
than jets as they do not require the observation of clearly structured “clusters” of parti¬ 
cles (jets). Contrary to jet production, a fragmentation function has to be introduced into 
the theoretical framework, cf. Eq. (Q), to take care of final-state collinear singularities. 
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FIGURE 4. As in Fig. ^but now for high-pj jet production. The “error bars”, Eq. (^, are for the STAR 
experiment taking into account its acceptance. Figure taken from [jlll]. 


In case of pion production, the are, however, fairly well constrained by data. It 
should be also emphasized that in the unpolarized case, the comparison between NLO 
theory predictions with jet production data from the Tevatron is extremely successful. 
The same is true for first preliminary data on the py.-spectrum for pions at \fs = 200 GeV 
from PHENIX [|g. 

The NLO QCD corrections to polarized jet production are available as a Monte Carlo 
code o- Apart from a significant reduction of the scale dependence, they are also 
mandatory for realistically matching the procedures used in experiment in order to group 
final-state particles into jets. For single-inclusive high-p^, hadron production the task of 
computing the NLO corrections has been completed only very recently [TSl] . Figure ^ 

shows for single-inclusive jet production at the NLO level as a function of the jet 
Pj- A cut in rapidity, |r]| < 1, has been applied in order to match the acceptance of 
STAR. The asymmetries turn out to be smaller than for prompt photon production, but 
thanks to the much higher statistics one can again easily distinguish between different 
spin-dependent gluon densities. Results for single-inclusive tjP production are presented 
in Fig. Note that here the expected statistical accuracy refers to only a very moderate 
integrated luminosity and beam polarization as targeted for the upcoming run of RHIC. 
Even under these assumptions a first determination of Ag can be achieved. 

The last process which exhibits a strong sensitivity to Ag is heavy flavor production. 
Here, the LO gluon-gluon fusion mechanism, gg —> QQ, dominates unless pj- becomes 
rather large. Unpolarized calculations have revealed that NLO QCD corrections are 
mandatory for a meaningful quantitative analysis. In the polarized case they have been 
computed recently in case of single-inclusive heavy quark production [jl^. Again, one 
observes a strongly reduced scale dependence for charm and bottom production at RHIC 
energies. It turns out that the major theoretical uncertainty stems from the unknown 
precise values for the heavy quark masses []T4|]. Since the heavy quark mass already 
sets a large scale, one can perform calculations for small transverse momenta or even 
for total cross sections which, in principle, give access to the gluon density at smaller 
jc-values than relevant for high-pj. jet, hadron or prompt photon production. 






FIGURE 5. As in Fig. |3] but now for high-pj. Tr'^-production. Note that here the statistical “error bars” 
5 All have been estimated by assuming only P = 0.4 and ^ = 7pb^^ in Eq. which is a realistic target 
for the next RHIC run. Figure taken from Jl^]. 


Heavy flavors are not observed direetly at RHIC but only through their deeay prod¬ 
uets. Possible signatures for eharm/bottom quarks at PHENIX are inelusive-muon or 
eleetron tags or /ie-coincidences. The latter provide a mueh better c/Z7-separation which 
is an experimental problem. In addition, lepton detection at PHENIX is limited to 
|y| < 0.35 and 1.2 < |y| < 2.4 for electrons and muons, respectively. Since heavy quark 
decays to leptons proceed through different channels and have multi-body kinematics, 
it is a non-trivial task to relate, e.g., experimentally observed py-distributions of decay 
muons to the calculated py.-spectrum of the produced heavy quark. One possibility is to 
model the decay with the help of standard event generators like PYTHIA [ |T5l ] by com¬ 
puting probabilities that a heavy quark with a certain (pj-^y) is actually seen within the 
PHENIX acceptance for a given decay mode. Eigure shows a prediction for the charm 
production asymmetry at PHENIX in NEO QCD for the inclusive-electron tag. The 
sensitivity to Ag is less pronounced than for the processes discussed above. It remains to 
be checked if heavy flavor production at RHIC can be used to extend the measurement 
of Ag towards smaller x-values. 


Further Information on Aq and Aq 

Inclusive DIS data only provide information on the sum of quarks and antiquarks for 
each flavor, i.e., Ag + Ag. At RHIC a separation of Am, Am, Ad, and Ad can be achieved 
by studying VT^-boson production. Exploiting the parity-violating properties of VT^- 
bosons, it is sufficient to measure a single spin asymmetry, , with only one of the 
colliding protons being longitudinally polarized. The idea is to study as a function 
of the rapidity of the VT, y^y, relative to the polarized proton [|T^. In EO it is then easy to 
show [11,1] that for VT+-production, ud —> VT+, and large and positive (negative) yy/,A^ 
is sensitive to Am/m (AJ/J). Similarly, -production probes Ad/d and Am/m. The NEO 












FIGURE 6. NLO single-inclusive charm production asymmetry (rescaled by 1 /jr™"') as a function of 
jc™" = for different sets of parton densities. The “error bars”, Eq. (^, are for the PHENIX 

experiment and include a detection efficiency for the channel c ^ eX. Eigure taken from 0 - 


QCD corrections for as well as the factorization scale dependence are small [[T^ . 
Experimental complications [Q] arise, however, from the fact that neither PHENIX nor 
STAR are hermetic, which considerably complicates the reconstruction of y^. Therefore 
it is important to understand on the decay-lepton level. Here, fully differential NEO 
cross sections are available as a MC code [j^ . The anticipated sensitivity of PHENIX 
on the flavor decomposed quark and antiquark densities is illustrated in Eig. ^ 

Semi-inclusive DIS measurements, ep —> HX, are another probe to separate quark and 
antiquark densities. HERMES has recently published first preliminary results [[T^]. The 
accessible x-range for the Aq and Aq densities is comparable to that of RHIC, see Eig. 
but at scales 2 ~ 1 — 2 GeV rather than Myy. The combination of both measurements can 
provide an important test of the QCD scale evolution for polarized parton densities and 
of the possible relevance of higher twist contributions at low scales. 


Towards a Global Analysis of Upcoming Data 

Having available at some point in the near future data on various different reactions, 
one needs to tackle the question of how to set up a “global QCD analysis” for spin- 
dependent parton densities. The strategy is in principle clear from the unpolarized case: 
an ansatz for the densities, Eq. (^, at some initial scale jUq, given in terms of some 
functional form with a set of free parameters, is evolved, Eq. (^, to a scale p relevant for 
a certain data point. A ;(;^-value is assigned that represents the quality of the comparison 
of the theoretical calculation to the experimental point. The parameters are varied until 
eventually a global minimum in is reached mutually for all data points. In practice, 
this approach is not fully viable since the numerical evaluations of the cross sections 
in NEO QCD are usually time-consuming as they require several tedious integrations. 
Hence the computing time for a QCD fit easily becomes excessive. 

In the unpolarized case, the wealth of DIS data already provides a pretty good knowl¬ 
edge of the parton densities, and reasonable approximations can be made for the most 
time-consuming processes. Eor instance, one can absorb all NEO corrections into some 
pre-calculated “correction factors” K, and simply multiply them in each step of the fit 
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FIGURE 7. Expected statistical accuracy for t^qjq from overlayed on two sets of parton densities. 
The full [open] circles refer to Aj^(lT+) [Aj^(VT^)]. Figure taken from 0. 


to the LO approximation for the eross seetions whieh ean be evaluated mueh faster. In 
the polarized ease, it is in general not at all elear whether sueh a strategy will work. 
Here, parton densities are known with much less aecuraey so far. It is therefore not pos¬ 
sible to use .fiT-faetors reliably. In addition, spin-dependent parton densities as well as 
partonie eross seetions may oseillate, i.e., have zeros, in the kinematieal regions of in¬ 
terest sueh that predietions at LO and the NLO ean show marked differenees. Clearly, in 
the polarized ease the goal must be to find a way of implementing efficiently, and with¬ 
out approximations, the exact NLO expressions for all relevant hadronic cross sections. 
A very simple and straightforward method based on “double Mellin transformations” 
was proposed in p^ ] . Recently, its actual practicability and usefulness in a global QCD 
analysis has been demonstrated pll ] in a case study based on fictitious prompt- 7 data. 


EXPLORING PHYSICS BEYOND THE STANDARD MODEL 


Spin observables are also an interesting tool to uncover important new physics. One idea 
is to study single spin asymmetries for large-pj. jets. In the standard model can 
be only non-zero for parity-violating interactions, i.e., QCD-electroweak interference 
contributions, which are fairly small. The existence of new parity-violating interactions 
could lead to sizable modifications Possible candidates are new quark-quark 


contact interactions, characterized by a compositeness scale A. RHIC is surprisingly 
sensitive to quark substructure at the 2 TeV scale, and is competitive with the Tevatron 


despite the much lower c.m.s. energy p2|] . Other candidates for new physics are possible 
new gauge bosons, e.g., a leptophobic Z'. Of course, high luminosity and precision 
as well as a good knowledge of polarized and unpolarized parton densities and of the 
standard model “background” are mandatory. For details, see 






SUMMARY AND OUTLOOK 


With first data from RHIC hopefully starting to roll in soon, we ean address many open, 
long-standing questions in spin physios like the longitudinally polarized gluon density. 
With data from many different prooesses taken at high energies, where perturbative QCD 
should be at work, a first global analysis of spin-dependent parton densities will be 
possible. For a long time to eome RHIC will provide the best souree of information on 
polarized parton densities, eertainly mueh improving our knowledge of the spin struoture 
of the nuoleon, and, perhaps, the next “spin surprise” is just round the oorner. Future 
projeets like the EIC []23|], whieh is ourrently under serutiny, would help to further deepen 
our understanding by probing aspeets of spin physios not aeoessible in hadron-hadron 
oollisions. The strueture funetion gj at small x or the spin eontent of oiroularly polarized 
photons are just two examples. 
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